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We have calculated the large-g expansion for the energy and magnetization cumulants at the first order phase 
transition point in the two-dimensional q-state Potts model to the 21st or 23rd order in 1/y/q using the finite 
lattice method. The obtained series allow us to give highly convergent estimates of the cumulants for q > 4. The 
results confirm us the correctness of the conjecture by Bhattacharya et al. on the asymptotic behavior of the 
energy cumulants for q — *• 4 + and a similar new conjecture on the magnetization cumulants. 



1. INTRODUCTION 

The g-state Potts model [[l]j|] in two dimen- 
sions has been investigated intensively as the test 
ground for analyzing the phase transition in many 
physical systems. Especially it is interesting be- 
cause the order of the phase transition changes 
when the parameter q is varied, i.e., from the 
first order for q > 4 to the second order for q < 4. 
The amplitudes of many quantities at the first or- 
der transition point are known exactly, including 
the latent heat, the spontaneous magnetization^ 
and the correlation length 0-|^). Bhattacharya 
et made a stimulating conjecture on the 

asymptotic behavior of the energy cumulants (in- 
cluding the specific heat) at the first order tran- 
sition point. Bhattacharya et also made the 
large-g expansion of the energy cumulants to or- 
der 10 in z = l/^/q and analyzing the expansion 
series with the conjecture they could give the es- 
timates of the cumulants at the transition point 
for q > 7 that are better than those given by 
other methods including the Monte Carlo sim- 
ulations^] and the low-(and high-)temperature 
expansions [p~0| p"2[ . The series obtained by Bhat- 
tacharya et al. are, however, not long enough to 
investigate the behavior of the energy cumulants 
for q closer to 4. 

Here using the finite lattice method |L3]-[l5) we 
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have generated the large-g series for the energy 
cumulants and the magnetization cumulants at 
the transition point to order 21 or 23 in z[^6|. The 
finite lattice method can in general give longer 
series than those obtained by the diagrammatic 
method especially in lower space (and time) di- 
mensions. In the diagrammatic method, one has 
to list up all the relevant diagrams and count the 
number they appear. In the finite lattice method 
we can skip this job and reduce the main work to 
the calculation of the partition function for a se- 
ries of finite size lattices, which can be done with- 
out the diagrammatic technique. This method 
has been used mainly to generate the low- and 
high-temperature series in statistical systems and 
the strong coupling series in lattice gauge theory. 
One of the purposes of our work is to demonstrate 
that this method is also applicable to the series 
expansion with respect to the parameter other 
than the temperature or the coupling constant. 

2. ENERGY CUMULANTS 

The latent heat C and the correlation length £ 
at the transition point are known to vanish and 
diverge, respectively, at q — > 4 + as 

C ~ 3vnr- 1/2 , (1) 



Table 1 

The specific heat for some values of q. The exact 
correlation length is also listed. 



q 


c d 


C 


u 


5 


2889(2) 


2886(3) 


2512.2 


6 


205.93(3) 


205.78(3) 


158.9 


7 


68.738(2) 


68.513(2) 


48.1 


8 


36.9335(3) 


36.6235(3) 


23.9 


9 


24.58761(8) 


24.20344(7) 


14.9 


10 


18.38543(2) 


17.93780(2) 


10.6 


12 


12.401336(3) 


11.852175(2) 


6.5 


15 


8.6540358(4) 


7.9964587(2) 


4.2 


20 


6.13215967(2) 


5.36076877(1) 


2.7 




Figure 1 . The behavior of the ratio of the second 

(2) 

energy cumulants to i 



with x = exp (-7T I Iff) and 2 cosh (9 = Jq. Bhat- 
tacharya et aVs conjecture says that the n-th en- 
ergy cumulants i^i" at the first order transition 
point P = Pt will diverge at q — ► 4 + as 

ff',(-l)7W ~ QB »-2£^_zii a; 3n/2-2 _ (3) 

This is from the fact that for /3 — ► fit at q = 4 
(the second order phase transition) the correla- 
tion length and the second cumulants diverge as 
£ ~ A|/3 - &|- 2 / 3 and F< 2 ) - j*|/3 - Al" 2/3 , re- 
spectively, so that - ^ r w2?^ (g/A) 3n/2 ~ 2 

and from the assumption that this relation be- 
tween the correlation length and the energy cu- 
mulants are also kept for q — > 4 + with (3 = j3 t . 
The constants a and -B in Eq.(||) should be com- 
mon to the ordered and disordered phases from 
the duality relation for each n-th cumulants. 

If this conjecture is true, the product 
p(n)£3n-4 j g & sm00 th function of 9, so we can 
expect that the Pade approximation of F^C P 
will give convergent result at p = 3n — 4. It 
has been tried for the large-g series obtained by 
the finite lattice methods for n = 2, • • • ,6 both 
in the ordered and disordered phases, which in 
fact gives quite convergent Pade approximants 
for p = 3n — 4 and as p leaves from this value 
the convergence of the approximants becomes bad 
rapidly. We give in Table 1 the values of the spe- 
cific heat C — PfF^ 2 ' evaluated from these Pade 



approximants for some values of q and present 
in Fig.l the behavior of the ratio F^/x plotted 
versus 9. These estimates are three or four orders 
of magnitude more precise than (and consistent 
with) the previous result for q > 7 from the large- 
q expansion to order z 10 by Bhattacharya et ai.Q 
and the result of the Monte Carlo simulations for 
q = 10, 15, 20 carefully done by Janke and Kap- 
plerQ. What should be emphasized is that we 
obtained the values of the specific heat in the ac- 
curacy of about 0.1 percent at q — 5 where the 
correlation length is as large as 2500. As for the 
asymptotic behavior of F^ at q — > 4 + , the Pade 

(2) (2) 

data of Fj fx and F D fx shown in Fig. 2 have the 
errors of a few percent around q = 4 and their be- 
haviors are enough to convince us that the conjec- 
ture (|) is true for n = 2 with a = 0.073 ± 0.002 . 
Furthermore from the conjecture (Jfy the com- 
bination {r(n-f) |F(")|/r(|)F( 2 )}^a;-i is 
expected to approach the constant B for each 
n(> 3), and in fact the Pade data for every 
n(= 3, • • • , 6) gives B — 0.38 ± 0.05 , which gives 
strong support to the conjecture also for n > 3. 

3. MAGNETIZATION CUMULANTS 

The behavior of the n-th magnetization cumu- 
lants M<") for /3 -► f3 t at q = 4 is well known 
as m!" ~ A? 1 ' (£) t™~ 2 and we can make a con- 
jecture parallel to the conjecture for the energy 
cumulants by Bhattacharya et ai: this relation 



Table 2 

The magnetic susceptibility for some values of q. 



q 


Xd 


Xo 


5 


9.13(3) x 10 4 


9.01(3) x 10 4 


6 


6.585(4) x 10 2 


6.665(4) x 10 2 


7 


70.54(1) 


77.31(1) 


8 


19.359(1) 


21.525(1) 


9 


8.0579(1) 


9.0106(2) 


10 


4.23276(2) 


4.73823(4) 


12 


1.720645(2) 


1.918274(3) 


15 


0.7304214(1) 


0.8056969(2) 


20 


0.309365682(1) 


0.33556421(1) 



cumulants and the correlation length for q = 4 
and (3 — >• fit (the second order phase transition) 
is kept in their asymptotic behavior for q — > 4 + 
at f3 = fit (the first order transition point). If 
this kind of relation is satisfied for the quantities 
at the first order phase transition point in more 
general systems as the asymptotic behavior when 
the parameter of the system is varied to make the 
system close to the second order phase transition 
point, it would serve as a good guide in investi- 
gating the property of the system. 
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